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Abstract: The unreasonable nitrogen (N) supply and low productivity are the main factors restricting the 
sustainable development of processing tomatoes. In addition, the mechanism by which the N application 
strategy affects root growth and nitrate distributions in processing tomatoes remains unclear. In this study, 
we applied four N application levels to a field (including 0 (NO), 200 (N200), 300 (N300), and 400 (N400) 
kg/hm/?) based on the critical N absorption ratio at each growth stage (planting stage to flowering stage: 
22%; fruit setting stage: 24%; red ripening stage: 45%; and maturity stage: 9%). The results indicated that 
N300 treatment significantly improved the aboveground dry matter (DM), yield, N uptake, and nitrogen 


use efficiency (NUE), while N400 treatment increased nitrate nitrogen (NO3-N) residue in the 20-60 cm 
soil layer. Temporal variations of total root dry weight (TRDW) and total root length (TRL) showed a 
single-peak curve. Overall, N300 treatment improved the secondary root parameter of TRDW, while 
N400 treatment improved the secondary root parameter of TRL. The grey correlation coefficients 
indicated that root dry weight density (RDWD) in the surface soil (0—20 cm) had the strongest relationship 
with yield, whereas root length density (RLD) in the middle soil (20—40 cm) had a strong relationship with 
yield. The path model indicated that N uptake is a crucial factor affecting aboveground DM, TRDW, and 
yield. The above results indicate that N application levels based on critical N absorption improve the 
production of processing tomatoes by regulating N uptake and root distribution. Furthermore, the results 
of this study provide a theoretical basis for precise N management. 
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1 Introduction 


Globally, the use of nitrogen (N) fertiliser has contributed to alleviating the problem of food 
security, particularly in realising high crop yields (Wang et al., 2018; Matiwane et al., 2019). 
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China is one of the countries with the largest N application in the world, the unreasonable 
application of N fertiliser is prevalent (Li and Lin, 2021). However, increasing the input of N 
fertiliser is unlikely to increase yield sustainably due to the decline in nitrogen use efficiency 
(NUE) (Liu et al., 2019; Liang et al., 2020). Moreover, blind N application leads to soil 
acidification and groundwater pollution in farmlands (Xiang et al., 2007; Gu et al., 2019). 
Therefore, it has become imperative for sustainable agricultural development to adopt scientific N 
fertiliser application practices, improve NUE, and reduce N fertiliser pollution in the environment 
(Li and Lin, 2021). 

Critical N absorption is an important parameter for evaluating NUE in certain growth periods 
(Ulrich, 1952). In studies conducted in pastures without N limitation, the relationship between 
critical N concentration and aboveground dry matter (DM) has been established. The critical N 
concentration was defined after continuous verification and parameter modification and has been 
applied to most C3 and C4 plants (Greenwood et al., 1990; Wang et al., 2013). Through N 
nutrition diagnosis, field application rates of 240-320 kg/hm? N have been shown to satisfy rice 
plant growth (Ata-Ul-Karim et al., 2014). Although the critical N concentration model provides 
the total N application rates, the critical N application rates at each growth stage have not been 
implemented in agricultural practice. Therefore, the key to accurate N application is to determine 
the N application rates at each growth stage based on the critical N concentration model. 

Processing tomato is among the vegetable crops with large cultivation areas globally (Farneselli 
et al., 2015) and contains important nutrients the human body requires (Lu et al., 2019). China has 
become a major producer and supplier of processing tomatoes, second only to the United States 
(Elia and Conversa, 2012). In recent years, drip irrigation technology has developed rapidly in 
China, significantly impacting the root architecture of processing tomatoes compared with 
traditional flood irrigation (Wang et al., 2022). Combined with the advantages of drip irrigation 
under mulch, the coupling mechanism of N absorption and the plant growth of processing tomatoes 
is vital for dryland agriculture development (Zotarelli et al., 2009; Li et al., 2016; Kang et al., 2017). 

The spatial distribution of plant roots plays a crucial role in soil-shoot—atmosphere interactions 
(Liu et al., 2018). Soil nutrients and moisture become heterogeneous due to the interaction 
between roots, which impacts the community structure and composition (Day et al., 2003). 
Conversely, the spatial distribution of roots responds differently to soil nutrient and moisture 
conditions (Sun et al., 2018). For example, studies have shown that soil moisture altered root 
architecture and biomass of soybean (Fenta et al., 2014), while N deficiency resulted in notable 
changes in lateral roots of cotton (Zhu et al., 2022). In this study, we assume that the N 
application strategy improves the productivity of processing tomatoes by regulating root growth 
and nitrate distribution. Under drip irrigation under mulch, critical N absorption of processing 
tomatoes was used to calculate the critical N demand at each growth stage (Tei et al., 2002). 
Therefore, this research aimed to (1) investigate the effects of N application levels on the 
productivity of processing tomatoes, N absorption, and nitrate nitrogen (NO3-N) residue; (2) 
explore the dynamic root distribution of processing tomatoes under different N application levels; 
and (3) clarify the regulatory mechanism of N application based on critical N absorption model on 
the productivity and root distribution of processing tomatoes. 


2 Materials and methods 


2.1 Experiment site 


The field experiment was conducted from April 2018 to July 2019 at the agricultural experimental 
station of Shihezi University (44°21'N, 86°04'E; 450 m) in Shihezi City, Xinjiang Uygur 
Autonomous Region of China. The precipitation was 76 and 114 mm during the growing season 
of processing tomatoes in 2018 and 2019, respectively, mainly concentrated in May. The average 
daily temperatures during growing season in 2018 and 2019 were 24.5°C and 23.1°C, respectively 
(Fig. 1). The chemical properties of topsoil (0-40 cm) were shown in Table 1. 
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Fig. 1 Daily precipitation, maximum temperature (Tmax), and minimum temperature (Tmin) during the growing 
season of processing tomatoes in 2018 (a) and 2019 (b) 


Table 1 Initial chemical properties of topsoil (0—40 cm) in 2018 and 2019 


Year pH Total N (g/kg) Available P (g/kg) Available K (g/kg) Organic matter (g/kg) 
2018 7.53 0.85 0.02 0.34 10.12 
2019 7.50 0.76 0.02 0.26 12.35 


Note: N, nitrogen; P, phosphorus; K, potassium. 
2.2 Experimental design and field management 


Wang et al. (2013) developed a statistical method for constructing the critical N concentration 
dilution curve of processing tomatoes based on DM in northern Xinjiang. The critical N 
concentration dilution curve was calculated as follows: 


Ne = 4.352DM °?”4, (1) 
where Nc is the N concentration (%). 
Critical N uptake was calculated by the following formula: 


Nupt = 43.521DM””° , (2) 


where Nupt is the N uptake (kg/hm?). The proportion of the N application rate at each growth 
stage to the total N application rate was determined using N uptake (Jing et al., 2020). We 
calculated the increase in N uptake at each growth stage according to the increase in DM of 
processing tomatoes at that specific stage, representing an increase in the critical N application 
rate. Therefore, the proportion of the N application rate at each growth stage to the total N 
application rate in the entire growth period was determined as follows: planting stage to flowering 
stage: 22%; fruit setting stage: 24%; red ripening stage: 45%; and maturity stage: 9%. 

The experiment followed a completely randomised block design with three replicates. We set 
four N application levels, including 0, 200, 300, and 400 kg/hm? (i.e., NO, N200, N300, and 
N400, respectively). We used N fertiliser as the top dressing based on the proportion for each 
growth stage. Additionally, 210 kg/hm? of phosphate fertiliser and 150 kg/hm? of potassium 
fertiliser were simultaneously applied to the soil as base fertilisers. 
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Before the experiment, healthy seeds of processing tomatoes (Rieger 87-5) were sown in pots 
and provided with sufficient water. The seedlings were transplanted to the field on 29 April 2018 
and 26 April 2019, when four true leaves appeared. The arrangement of the drip lines was set as 
‘one mulch, two drip lines, and two rows of processing tomatoes’, with a plant spacing of 30 cm 
and a row spacing of 60 cm (Fig. 2). The area of each experimental plot was 86.4 m? (7.2 mx12.0 
m). Irrigation was performed every 7 to 10 d, and each irrigation was controlled using a water 
meter. Irrigation was postponed in the event of rainfall. The total irrigation volume during the 
entire growth period was approximately 4700 m*/hm/. 
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Fig. 2 (a), overhead view of the planting pattern of processing tomatoes in the experimental field; (b), 
stereoscopic view of root samples taken using the layered-mining method 


2.3 Sampling and measurements 


2.3.1 Measurement of aboveground dry matter (DM), nitrogen (N) uptake, yield, and nitrogen use 
efficiency (NUE) 
At the harvest of processing tomatoes in 2018 and 2019, three representative plants were randomly 
collected from each plot to determine aboveground DM and N uptake. The aboveground parts of 
the plants were divided and oven-dried at 75°C to a constant weight, and aboveground DM was 
determined through gravimetry. Subsequently, the oven-dried plant samples were pulverised and 
screened through a 0.1 mm sieve. An appropriate amount of fine powder was digested with 
H2S04-H20: to determine N concentration using the Kjeldahl method. N uptake was calculated as 
the product of DM and N concentration. The processing tomatoes were manually harvested from a 
2.4 m? area per plot to determine yield. NUE (%) was calculated as follows: 

NUE = NO «100%, (3) 

Na 

where Cy is the N uptake by plant in the areas where N is applied (kg/hm7); Co is the N uptake by 
plant in the areas without N application (kg/hm”); and Na is the application amount of N fertiliser 
(kg/hm’). 
2.3.2 Soil nitrate nitrogen (NO; -N) content 
After the harvest of processing tomatoes, soil samples were collected from 0-80 cm layers at 20 
cm intervals using an earth drill with a 2.5 cm diameter. The soil samples were dried at room 
temperature, and then milled and screened. The NO3 -N content was determined using a 1 mol/L 
KCI solution through the ultraviolet spectrophotometer method. 


2.3.3 Root sampling and analysis 


In each plot, the roots of processing tomatoes were regularly harvested using the excavation 
method. Because the roots were mainly distributed in the 0—60 cm soil layer, and 30 cm around 
the plant was the concentration area of water and nutrients under drip irrigation (Jing et al., 2021), 
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the roots were collected in a soil column (30 cmx30 cmx60 cm) at 20 cm intervals (Fig. 2b). The 
root samples were obtained using the manual flushing method and scanned with a special root 
scanner (Epson Perfection V700 Photo, EPSON, Nagano, Japan). The root length was determined 
by analysing the scanned images using the Winrhizo Pro Vision 5.0 analysis program. Finally, the 
root samples were dried at 75°C for 12 h to obtain the root dry weight. 

Root dry weight and root length per soil layer in the unit area were expressed by the following 
formulas: 


D 

RDW=—, 4 
5 (4) 
L 

RL=—, 5 
5 (5) 


where RDW is the root dry weight (g/m”); D is the root dry weight per soil layer (g); S is the soil 
area (m°); RL is the root length (m/m7”); and L is the root length per soil layer (m). Total root dry 
weight (TRDW; g/m) or total root length (TRL; m/m?) was the sum of root dry weight or root 
length of each soil layer. 

Similarly, root dry weight density as well as root length density in unit volume of soil were 
expressed by the following formulas: 


RDWD=2., (6) 
V 
L 

RLD=—_, 7 
7 (7) 


where RDWD is the root dry weight density (g/m*); RLD is the root length density (m/m°*); and V 
is the soil volume (m°). 
2.3.4 Quantitative analysis of beta model 


The beta function of the root growth rate can determine the relationship between growth rate and 
duration (Yin et al., 2002; Wang et al., 2014). The formula is as follow: 


te 

= KEAMEN A 28 ET, 8 

W = Wmax i USEK 21, —bys Y< mS le» (8) 
l= b t 


e 
where w is the growth parameters of root (TRDW or TRL); ¢ is the number of days after planting 
(d); Wmax is the maximum growth parameters of root when it grows to te day (g/m? or m/m? for 
TRDW or TRL, respectively); and tn is the number of days when the root growth rate reaches the 
maximum value (d). 

The average rate of root growth (ca; (g/m°)-d™! or (m/m?°)-d™! for TRDW or TRL, respectively) 
was expressed by the following formula: 


Ca = (9) 


The maximum root growth rate (Cm; (g/m7)-d! or (m/m?).d! for TRDW or TRL, respectively) 
was expressed by the following formula (Yin et al., 2002): 


bm 


2t, -t ty Metn 
C, =— "x G ) X Wnax - (10) 
— t 
te (te tm) 


2.3.5 Gray correlation 
The gray correlation was evaluated by the method of Jia et al (2020): 
Step 1: the original data transformation by equalization: 
X; 
a (11) 


oe 
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where x'; is the equalized transformation value; x; is the actual value of index; and X; the 


ij y 
average value of index j. 

Step 2: calculation of the correlation coefficient. The yield of processing tomatoes was used as 
the basic sequence, and the root growth parameters in each soil layer were set as the subsequence. 
The correlation coefficient was calculated by the following equation: 

Ly(K) = Amin + PAmax , (12) 
Ag; (K) + pA 
where Loi(K) is the correlation coefficient; Amin and Amax are the minimum and maximum values 
of the absolute differences at any moment, respectively; p is the resolution ratio, with the value 
set as 0.5; and Ao(K) is the absolute difference between the two comparison sequences at moment 
K. 
Step 3: correlation calculation. We used the following equation: 


1 N 
hi y rK), (13) 


where ro; is the correlation between the subsequence i and the basic sequence; and N is the length 
of the comparison sequences. 

Step 4: testing the significance of the correlation coefficient. The t value was calculated using 
the following equation: 


max 


t=—ixJn-2, (14) 


l-r; 
where n-2 is the degree of freedom. 
2.4 Statistical analyses 


The analysis of variance was conducted using SPSS 17.0 software. Path analysis was employed to 
assess the relationships between indicators using SPSS 17.0 software. The least significant 
difference tests were used for significance analysis. Plotting was completed using Origin 9.0 
software. 


3 Results 


3.1 Aboveground DM, yield, N uptake, and NUE under different N application levels 


The N application levels significantly affected the aboveground DM, yield, N uptake, and NUE of 
processing tomatoes (P<0.01; Fig. 3). Although there were differences in the results between the 
two years due to climatic factors, the same trends were observed in different treatments for both 
years. The order of aboveground DM, yield, and N uptake for the four N application levels was 
ranked as NO<N200<N400<N300. Furthermore, NUE was higher under N300 treatment 
compared to N200 and N400 treatments. 


3.2 Dynamic changes in total root dry weight (TRDW) and total root length (TRL) 


The TRDW and TRL showed similar temporal variations in the two years (Fig. 4). The growth 
parameters of roots were fitted using the beta model, with the determination coefficients (R°) 
exceeding 0.85, indicating that the beta model could effectively reflect root development. The 
growth process of roots under different N application levels showed a single-peak curve, 
characterised by slow growth, rapid growth, and then a gradual decline. There were differences in 
the TRDW and TRL under different N application levels (Table 2). During the two years, the 
order of the maximum growth parameters of root, the average rate of root growth, and the 
maximum root growth rate in TRDW for the four N application levels was 
NO<N200<N400<N300. Additionally, these values in TRL under N400 treatment were higher 
compared to other treatments. 
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Fig. 3 Aboveground dry matter (DM; a), yield (b), nitrogen (N) uptake (c), and nitrogen use efficiency (NUE; d) 
under different N application levels in 2018 and 2019. NO, N200, N300, and N400 represent the N application 
levels of 0, 200, 300, and 400 kg/hm?, respectively. Different lowercase letters in the same year indicate 
significant differences among different N application levels at P<0.05 level. ** indicates significance at P<0.01 
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Fig. 4 Beta model of the variations in total root dry weight (TRDW; a and b) and total root length (TRL; c and 
d) under different N application levels in 2018 and 2019. The fitting curves are the beta equation fitting curves of 


the N application levels at 0, 200, 300, and 400 kg/hm?. 
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Table 2 Secondary root parameters under different N application levels 


Secondary root parameter 
Root Ty p 


Year parameter N@PPHeation level — wa(gm? gy gy CaCl dT or cn (Ed or 
or m/m’)* e m (m/m’°)-d 1y (m/m?’)-d 
NO 134.344 86.62° 53.80° 1.55¢ 2.59% 
N200 216.16° 84.99” 52.25° 2.54° 4.21° 
TRDW 
N300 403.28" 81.60° 60.47* 4.948 10.19° 
N400 267.95° 79.90° 56.45° 3.35? 6.40° 
2018 
NO 355.76° 69.16° 43.05° 5.14° 8.59° 
TRE N200 408.77” 68.93° 41.14” 5.91° 9.61° 
N300 581.57" 71.61° 42.61” 8.12° 13.14° 
N400 596.53° 72.02? 47.24" 8.28° 14.48° 
NO 112.074 88.03* 47.97° 1.274 1.974 
N200 177.91° 81.47° 48.42° 2.02°¢ 3.535 
TRDW 
N300 308.10° 79.61” 57.08° 3.87° T39 
RT N400 222.65° 80.83° 56.17* 2.75° 5.14” 
01 
NO 348.194 71.68° 47.20° 5.36” 8.53” 
N200 371.98° 71.23* 45.19° 5:22? 8.86° 
TRL 
N300 514.03” 73.25° 42.58” 7.02* 11.20° 
N400 555.99° 74.85° 47.21* 7.43° 12.54° 


Note: NO, N200, N300, and N400 represent the N application levels of 0, 200, 300, and 400 kg/hm?, respectively. TRDW, total root dry 
weight; TRL, total root length; Wmax, the maximum growth parameters of root; te, the number of days when the root growth parameter 
reaches the maximum value; fm, the number of days when the root growth rate reaches the maximum value; ca, the average rate of root 
growth; Cm, the maximum root growth rate. Different lowercase letters in the same year indicate significant differences among different 
N application levels at P<0.05 level.” represents that the unit of Wmax is g/m’ for TRDW and m/m? for TRL; * represents that the 
unit of c, and Cm is (g/m?)-d"! for TRDW and (m/m?)-d"! for TRL. 

3.3 Soil NO3-N content, root dry weight density (RDWD), and root length density (RLD) 


in soil profiles 


At the harvest stage, the soil NO3-N content showed an initial increase followed by a decrease 
with an increase in soil depth under all treatments (Fig. 5). The order of the total soil NO3 -N 
content for the four N application levels was NO<N200<N300<N400. Remarkably, N400 
treatment resulted in a significant accumulation of soil NO3 -N content in the 20—60 cm soil layer. 
In addition, the vertical distribution of the roots showed adaptive changes. Both RDWD and RLD 
decreased gradually with an increase in soil depth (Fig. 6). The root was mainly distributed in the 
0-20 cm soil layer. Within the 0-20 cm soil layer, the order of RDWD and RLD under the four N 
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Fig. 5 Soil nitrate nitrogen (NO3-N) content in the 0-80 cm soil layers under different N application levels in 
2018 (a) and 2019 (b). Bars are standard errors. 
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Fig. 6 Changes in root dry weight density (RDWD; a and b) and root length density (RLD; c and d) in different 
soil layers under different N application levels in 2018 and 2019. Bars are standard errors. 


application levels was NO<N200<N400<N300 and NO<N200<N300<N400, respectively. In the 
20—40 and 40-60 cm soil layers, the RDWD and RLD under N300 treatment were higher 
compared to other treatments. 


3.4 Grey correlation coefficients between yield and root density in different soil layers 


Under N application levels based on critical N absorption, the minimum grey correlation 
coefficient between yield and RDWD or RLD was observed in the 40—60 cm soil layer. The 
maximum grey correlation coefficient between yield and RDWD was found in the 0-20 cm soil 
layer, while the maximum grey correlation coefficient between yield and RLD was observed in 
the 20—40 cm soil layer (Table 3). Additionally, N application levels affected the grey correlation 
coefficient in the 0-60 cm soil layer. For RDWD, the grey correlation coefficient of N300 
treatment was higher than that of other treatments in different soil layers. For RLD, the grey 
correlation coefficient of N400 treatment was higher than that of other treatments in the 0-20 and 
40-60 cm soil layers, while the grey correlation coefficient of NO treatment was higher than that 
of other treatments in the 20—40 cm soil layer. These results indicated that RDWD in the surface 
soil (0-20 cm) had the strongest relationship with yield, but its influence decreased with 
increasing soil depth. The RLD in the 20-40 cm soil layer also had a strong relationship with 
yield, significantly impacting yield under N deficiency (P<0.05). 


3.5 Path analysis of aboveground DM, TRDW, N uptake, soil NO3 -N content, and yield 


The path analysis model was used to assess the relationships among aboveground DM, TRDW, N 
uptake, soil NO3 -N content, and yield (Fig. 7). The direct impact degree on yield followed the 
order of N uptake>TRDW>aboveground DM>soil NO3-N content. Moreover, the coefficients of 
N uptake and TRDW on yield were significant (P<0.01). N uptake significantly and positively 
affected TRDW and aboveground DM (P<0.05). In addition, aboveground DM and soil NO3 -N 
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content were negatively related to yield. These findings indicated that an increase in nutrient 
growth and soil N does not necessarily translate into continuous yield improvement and may 
result in plant overgrowth and N leaching. 


Table 3 Grey correlation coefficients between root parameters in different soil layers with yield under different 
N application levels 


Soil depth (cm) 


N application 


Root parameter level 2018 2019 
0-20 20-40 40-60 0-20 20-40 40-60 
NO 0.649* 0.536 0.494 0.614 0.548 0.599 
naw N200 0.651" 0.501 0.515 0.729" 0.544 0.587 
(g/m’) N300 0.685 0.645 0.595 0.759" 0.677 0.657 
N400 0.635 0.534 0.478 0.693 0.549 0.523 
NO 0.559 0.859" 0.557 0.641 0.727* 0.554 
RED N200 0.573 0.676 0.508 0.583 0.637 0.540 
(m/m’) N300 0.553 0.559 0.581 0.499 0.576 0.481 
N400 0.673 0.757" 0.610 0.671 0.726 0.663 


Note: RDWD and RLD are root dry weight density and root length density, respectively. * and ™ indicate significant correlation 
coefficient at P<0.05 and P<0.01 levels, respectively. 


0.911" 


Aboveground DM = 


Soil NOs-N content 


0.811" a 0.563"* 


Fig. 7 Path model among aboveground DM, TRDW, N uptake, soil NO3-N content, and yield. Blue lines 
represent positive relationships, while red lines represent negative relationships. The width of the arrows 
represents the strength of significant standardised path coefficients. * and “ indicate significant path coefficient at 
P<0.05 and P<0.01 levels, respectively. 


4 Discussion 


4.1 Productivity and N uptake responses to various N application levels 


Reasonable fertilisation plans have been proven effective in promoting aboveground DM and 
yield (Ronga et al., 2017; Banger et al., 2020). In this study, we implemented an N application 
strategy based on critical N absorption at each growth stage. Our results showed a positive 
correlation between N application rate and aboveground DM or yield of processing tomatoes 
within 300 kg N/hm*. However, the 400 kg/hm? N application hindered aboveground DM and 
yield, consistent with the findings of Farneselli et al. (2018). This could be caused by the 
following reasons. First, an adequate N input ensures root growth and water absorption, thereby 
promoting aboveground DM and yield. However, excessive N input causes inorganic N 
accumulation in the soil, resulting in soil salinisation (Cheng et al., 2021). Second, high N 
concentration directly reduced soil microbial activity by enhancing specific ion toxicity, thereby 
inhibiting root growth and nutrient absorption (Yan et al., 2021). 

Reasonable regulation of soil N holds considerable potential for maintaining stable yield and 
NUE in arid areas (Zheng et al., 2021). In this study, N300 treatment had the highest N uptake 
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and NUE, while N400 treatment had the highest soil NO3-N content. Despite the plastic film 
mulching increasing soil water storage, the high N concentration accelerated soil moisture 
reduction due to excessive vegetative growth and an increased transpiration rate (Long et al., 
2021). Therefore, under water-deficit conditions, the migration rate of NO3-N to deep soil could 
be slower in N400 treatment. In addition, high N concentration inhibits plant N uptake and 
increases NO3 -N residue (Yuan et al., 2018). 


4.2 Root growth responses to various N application levels 


In this study, the secondary root parameters were analysed using the beta model. Overall, 
compared with NO and N200 treatments, N300 treatment improved the secondary root parameter 
of TRDW, while N400 treatment improved the secondary root parameter of TRL. Wang et al. 
(2014) reported that the community structures of crops were optimised by maintaining high 
TRDW in the late growth period under a reasonable planting mode. These improvements may be 
due to nutrient deficiency under NO and N200 treatments, which led to decreased soil microbial 
activity and hindered root growth (Isobe et al., 2018). In addition, nutrient deficiency in the late 
growth period led to the blockage of plant physiological metabolism (Zaman et al., 2021). Most 
roots of processing tomatoes were accumulated in the uppermost (0—20 cm) soil layer. Generally, 
plant root configuration is controlled by genetic factors (Clark et al., 2011; Ning et al., 2019). 
However, plant roots exhibit adaptive changes in response to specific environmental conditions in 
highly distinct growth environments (Martins et al., 2019). 

Appropriate N concentration promoted root growth (Hackett, 1972; Drew, 1975; Postma et al., 
2014), while high N concentration inhibited it (Forde and Lorenzo, 2001). This phenomenon was 
also observed in different soil layers in this study. Specifically, the RDWD of N300 treatment was 
highest compared to other treatments in the 0—20 cm soil layer, while the RLD of N400 treatment 
was the largest. This could be attributed to the following two mechanisms. First, high N 
concentration inhibits the biosynthesis and transport of indole acetic acid, thereby inhibiting 
taproot growth (Walch-liu et al., 2001). Second, lateral root growth is hindered by the inhibition 
of dividing tissue activity (Zhang and Forde, 1998). Additionally, research has shown that high N 
concentration promotes the regeneration of superficial capillary roots, which may explain the 
longer root length in the topsoil of N400 treatment. In this study, the RLD of N300 treatment was 
highest compared to other treatments in the 40—60 cm soil layer. Remarkably, the increased RLD 
in the deeper soil layer effectively facilitated nutrient and water absorption (Zhang et al., 2009; 
Jia et al., 2020). 


4.3 Crucial determinants of promotion in productivity 


N is essential for normal plant metabolism and is the main limiting factor for crop yields (Shi et 
al., 2022). The distribution of carbohydrates among various plant organs ensures a balanced 
growth ratio between underground and aboveground parts, but applying N fertilisers can disrupt 
this growth balance (Lovelli et al., 2012). Thus, N fertilisers significantly affect the biomass of 
plant shoots and roots (Sainju et al., 2017). Furthermore, N fertilisation plays a crucial role in soil 
pore structure and microbial balance (Liang and Shi, 2021). Therefore, employing the path model 
is necessary to understand the relationships among N uptake, aboveground DM, TRDW, soil 
NO3-N content, and yield. The present path model indicated that N uptake is a crucial factor 
affecting aboveground DM, TRDW, and yield. Moreover, aboveground DM and soil NO3 -N 
content were negatively related to yield, suggesting increased vegetative growth and the risk of 
NO3-N leaching. Notably, no significant negative correlation was observed between 
aboveground DM and yield, mainly due to the vigorous nutritional growth of plants under N400 
treatment. However, N300 treatment obtained the highest yield. According to the nutrient status 
of the soil and plants, plants use specialised vascular tissue to transport N, effectively regulating 
biomass distribution (Hermans et al., 2006; Litton et al., 2007; McCarthy and Enquist, 2007). 
However, due to the soil water infiltration caused by rainfall and irrigation, N leaching was 
inevitable (Wang and Li, 2019). In conclusion, we confirmed the crucial role of N application 
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based on critical N absorption in improving productivity, although the risk of NO3-N leaching 
still exists. Moreover, further validation is needed to understand the physiological mechanisms 
underlying the responses of processing tomatoes to critical N absorption. 


5 Conclusions 


Under N application levels based on critical N absorption, N300 treatment significantly improved 
the productivity and NUE of processing tomatoes, while N400 treatment increased the risk of 
NO3 -N leaching. The roots under N300 and N400 treatments showed increased growth, which 
promoted yield. While the RDWD in the surface soil (0-20 cm) had the strongest relationship 
with yield, the RLD in the middle soil (20-40 cm) also had a strong relationship with yield. 
Additionally, N uptake was a crucial factor influencing aboveground DM, TRDW, and yield. Our 
results suggested that the 300 kg/hm? N application, based on the critical N absorption ratio 
(planting stage to flowering stage: 22%; fruit setting stage: 24%; red ripening stage: 45%; and 
maturity stage: 9%), improved root distribution and soil NO3-N content and promoted the 
productivity processing tomatoes. 
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